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Abstract: This paper describes a state estimation approach for non-causal time- varying linear 
descriptor equations with uncertain parameters. The uncertainty in the state equation and in 
the measurements is supposed to admit a set-membership description. The approach is based on 
the notion of the linear minimax estimation. Suboptimal minimax state estimation algorithm is 
introduced for DAEs with non-stationary rectangular matrices. Optimal algorithm is presented 
for DAEs with special structure of the matrices. A comparison of suboptimal and optimal 
algorithms is presented for 2D time- varying DAE with a singular matrix pencil. 

Keywords: Set-membership state estimation, descriptor systems, singular systems, DAEs, 
minimax. 



1. INTRODUCTION 

Dynamical systems described by coupled differential and 
algebraic equations (DAEs) arise naturally in many appli- 
cations. In particular, DAEs occur in econometrics (singu- 
lar dynamic Leotief systems Luenberger and A. (1977)), 
modelling of constrained multibody systems Mills (2006), 
electrical circuit synthesis Rcis (2008), bioprocess and 
chemical engineering Becerra et al. (2001), representing 
a common modelling tool Mehrmann and Stykel (2005). 
On can divide contributions to the theory 1 of DAEs 
into results for casual DAEs and results for non-causal 
systems. 1.1 Causal DAEs. The solvability theory for 
finite-dimcntional systems with constant coefficients 

Fx t = Cx + Bf (1) 

is based on the reduction of the matrix pencil sF — C to 
the Kroncckcr canonical form: if dct[sF — C] ^ then 
for all initial values x(to) = x there exists a unique 
solution x(-). Changing the basis in the state space and 
differentiating exactly s times (d is an index of the pencil 
sF — C), one can reduce (1) to some equivalent Ordinary 
Differential Equation (ODE), provided /(•) is sufficiently 
smooth and meets some algebraic constraints. The details 
of the reduction process are presented in Gantmacher 
(1960). The index d of the pencil sF — C is said Campbell 
and Pctzold (1983) to be an index of linear DAE (1). 

The notion of a Standard Canonical Form (SCF) allows 
to generalize the index approach on variable coefficients: 
for instance, in Campbell and Petzold (1983) it was shown 
that (1) with analitical F(-), C(-) and B(-) is solvable 2 iff 



1 Solvability and numerical methods, control and observation for 
DAE and so on. 

2 DAE (1) is said to be solvable at [io,T] if for every sufficiently 
smooth /(•) there exists at least one continuously differentiablc 
solution, provided F(-), C(-) are sufficiently smooth. 



(1) can be converted into SCF. In Campbell (1987) it was 
noted that not all solvable DAEs can be put into SCF and 
the solvable DAE is equal to some differential-algebraic 
equation in the canonical form which generalize SCF. In 
this regard, we say that the DAE is causal if 1) it can be 
reduced (at least locally in the non-linear case) to normal 
ODE and 2) if it is solvable for the given intial condition 
Xo and input /(•) then the solution is unique. 

The geometry of the reduction procedure for nonlinear 
causal DAEs F(x, x) = was investigeted in Reich (1990); 
Rabier and Rheinboldt (1994), where the index of DAE 
was defined as a smallest natural d so that the sequence 
of the constraint manifolds Reich (1990) 

M k := lW fe _i n M fc _i, W k = {x e K" : {x,p) E M k } 

with M := {(x,p) : F(x,p) = 0} becomes stationary 
for k > d. This coincides with the definition of the 
index of linear DAE. Solvability of the linear casual DAEs 
with impulses in the input was addressed in Rabier and 
Rheinboldt (1996). Further discussion of the linear DAEs 
solvability theory in finite-dimensions and related topics 
are presented in Mehrmann and Stykel (2005); Samoilcnko 
et al. (2000). 

Basic ideas of the index approach (reduction of the pencil 
sF — C to the canonical form) were extended on systems 
with constant operator coefficients in Rutkas (1975), pro- 
vided the poles of the operator-valued function s n- (sF — 
C) _1 are contained in some bounded vicinity of 0. 

1.2 Non-causal DAEs. The non-causal DAE may have 
several solutions. For instance, consider 

|[l 0](£)= [ci <*](£) + /(*),[! 0](Z)(t )=X 

Let x 2 (-) € L 2 (t ,T), /(•) G U(t ,T) and x Q E K. By 
definition put 



xi(t) := exp(ci(t-t ))a;o+ / exp(ci(t-s))c 2 X2(s)+ f (s)ds 

It is clear that any solution of this DAE is given by the 
formula t i->- (xi(t), x 2 (t)) T . According to a behavioral 
approach Ilchmann and Mehrman (2005) one can think 
about x 2 as an input or as a part of the system state. 
Reduction of the non-causal DAEs with variable matrices 
was studied in Shlapak (1975) and Eremcnko (1980), 
provided F(-) has a constant rank. In this case DAE may 
be splitted into differential and algebraic equation. Further 
splitting is possible under some restrictions on matrices of 
resulting system. 

In this paper we focus on the inverse problem for linear 
non-causal DAEs with rectangular non-stationary matri- 
ces: given measurements y(t), t G [io>^] of the DAE' 
solution x(-), to reconstruct x(-). Here x(-) is said to be a 
solution of 

(Fx) t (t) = C(t)x(t) + f(t) (2) 

with initial condition Fx(to) = 0, if Fx(-) is absolutely 
continuous function, (Fx)t belong to h 2 (to, T, R m ), x(-) 
verifies (2) almost everywere and Fx(to) = holds. 
This definition guarantees that a linear mapping D, in- 
duced by (2), is closed as a mapping from h 2 (t n ,T,R n ) 
to h 2 {t ,T,R m ) Zhuk (2007). This, in turn, allows to 
properly define the system Zhuk (2007), adjoint to (2), 
that is of primary interest in optimal control. In addition, 
the method of matrix pensils is sometimes difficult to 
apply in the finite-dimensional optimal control theory For 
instance Ozcaldiran and Lewis (1989), applying the linear 
proportional feedback / = Kx to (1) with regular pencil 
det(sF-C) one may arrive to the system with singular 
pensil det(sF — C — BK) = 0. Therefore it is reasonable 
to apply the above definition of the DAEs' solution in the 
control framework. A much general one is presented in 
Kurina and Marz (2007), where a properly stated leading 

term A(t)—F(t)x is used in order to give a feedback 
solution to LQ-control problem with DAE constraints. 

In what follows we assume that y(t) = H(t)x(t)+n(t), x(-) 
is a solution of (2) in the sense of the above definition, the 
noise ry(-) is a realization of a random process VP, the input 
/(•) is uncertain and belongs to the given set G. The aim is 
to construct a worst-case estimation of the inner product 
(£,Fx(T)), £ e R n as a function of y(-), assuming that G 
is bounded set and the correlation function of 'J belongs 
to the given bounded set R of matrix-valued functions. 
This problem was solved in Nakonechny (1978), provided 
F = I. The case of deterministic measurement's noise was 
addressed in Bcrtsckas and Rhodes (1971); Milanese and 
Tempo (1985); Kurzhanski and Valyi (1997), where the 
optimal worst-case estimation is shown to be a dynamical 
system, describing the evolution of the central point of the 
ODEs reachability set, consistent with observations. 

In this paper, we generalize the theory of minimax state 
estimation Nakonechny (1978) on a class of linear non- 
causal DAEs: F e R mxn and t ^ C(t) e R mxn C(t) 
is continuous on [io,? 1 ]. The same results can be proved 
for the deterministic noise r)(-) and ellipsoidal bound for 
uncertain /(•) and ??(•), giving the generalization of Bert- 
sekas and Rhodes (1971); Kurzhanski and Valyi (1997). 
The major contributions of this paper is an implemen- 



tation of the abstract Generalized Kalman Duality prin- 
ciple Zhuk (2009b) for non-causal time-dependent DAEs 
(Theorem 3). Duality allows to find and exact expression 
for the worst-case estimation error and to establish the 
necessary and sufficient conditions on I for the worst-case 
error to be finite. These conditions, in turn, defines some 
subspace C{T) in the state space of (2), which is called 
a minimax observable subspace. In fact, £(T) describes 
an "observable" (in the minimax sense) part of x{T) with 
respect to the measured y(t), t <G [to,? 1 ]: if ^ € C{T) then 
we can provide the worst-case estimation of (£, Fx(T)) 
with finite worst-case error (which describes the measure 
of how poor the estimation quality may be); otherwise 
the state x(T) is not observable in the direction £, that 
is for any estimation of (£,Fx(T)) the estimation error 
varies in [0, +oo], so that, for any linear estimation and 
natural N there is a realization of uncertain parameters 
/(•) and rj(-) such that the estimation error will be greater 
than N. Note, that the notion of the minimax observable 
subspace C(T) is a implementation (in the case of DAEs) 
of the abstract minimax observability concept, presented 
in Zhuk (2009b). Some aspects of classical observability for 
DAEs were considered in Campbell et al. (1991) for causal 
systems and in Frankowska (1990) for non-causal systems. 

As a result of application of Generalized Kalman Duality 
for non-causal non-stationary DAE we derive a suboptimal 
worst-case state estimation algorithm (Corollary 6). The 
algorithm gives a suboptimal estimation of the projection 
of the state x{T) onto the minimax observable subspace 
£(T). It is sequential, that is the algorithm is represented 
in terms of the unique solution to a Cauchy problem for 
some ODE, which has a realization of the observations 
y(t), t G [to,T] as the input. Therefore, it is sufficient to 
know measurements y(t), t e [T, Xi] and the estimation 
at t = T in order to compute the estimation of x(T\). 
The algorithm works for "non-Gaussian noise" rj unlike the 
family of Kalman-like estimators. The optimal algorithm 
is also presented (Corollary 9), provided the matrices of 
DAEs have "some regularity" (Proposition 7). Kalman 
filtering approach was previously applied to linear DAEs 
with constant coefficients in Gerdin et al. (2007); Darouach 
et al. (1997), provided sF — C is regular. In this regard 
we note that the latter assumption can be substituted 
by the less restrictive one: sF' — C — H'H is regular 
(see Example above). Further information on Kalman 
filtering for causal DAEs is to be found at Xu and 
Lam (2006). A worst-case state estimation for non-causal 
linear continuous DAEs with non-stationary rectangular 
matrices was not considered in the literature before. The 
notion of the minimax observable subspace was applied 
in Zhuk (2009a) in order to construct the optimal state 
estimation algorithm for discrete time non-causal DAEs. 

Notation: Ery denotes the mean of the random element rj, 
intG denotes the interior of G, /(•) or / denotes some 
element of the functional space, f(t) denotes the value 
of the function / at time t, h 2 (to , T, R m ) denotes the 
space of square-integrable functions with values in R m , 
H 1 (t , T, R m ) denotes the space of absolutely continuous 
functions with L 2 -derivative and values in R m , the super- 
script ' denotes the operation of taking an adjoint, c(G, ■) 
denotes the support function of some set G, S(G, •) denotes 



the indicator 3 of G, dom/ = {x : f(x) < oo}; (•,•) 
denotes the inner product in Hilbert space H, K™ denotes 
n-dimensional Euclidean space over real field, S > means 
(Sx, x) > for all x from within appropriate Hilbert space, 
L' denotes adjoint operator, F' denotes transposed matrix, 
F + denotes pseudoinverse matrix. 

2. LINEAR MINIM AX ESTIMATION FOR DAES 



A function u(y) = J t (u(t),y(t))dt is called an a priori 
minimax mean-squared estimation in the direction I (£- 
estimation) if inf„ a(T, u, £) = a(T,u,£). The number 
<t(T, £) = o~(ii, £) is called a minimax mean-squared a priori 
error in the direction £ at time- instant T (^-error). The 
set C{T) = {£ E K™ : &(£) < +00} is called a minimax 
observable subspace. 



Consider a pair of systems 

(Fx) t (t) = C(£)x(£) + f(t),Fx(t ) = 0, 

y(t) = H(t)x(t)+ V (t),te[t ,T], {6 > 

where x(t) G R™, /(£) G R m , y{t) G R p , rj(t) G W 
represent the state, input, observation and observation's 
noise respectively, F G W nxn , /(•) G L 2 (t ,T), C(t) and 
H(t) are continuous matrix-valued functions, to,F G M. 

According to Zhuk (2007) we say that x(-) is a solution 
of (3) if Fx(-) G H^io, T,M m ) and the derivative (Fx) t 
of Fx(-) coincides with the right side of (3) almost every- 
where (a.e.) on [to,T] and Fx (to) = 0. 

i?em«d 1. As F e R mxn the pencil F - AG(t) is ir- 
regular Muller (1998) implying (3) may have non-unique 
solution. In what follows we will refer such DAEs as non- 
causal Muller (1998). 

In the sequel we assume rj(-) is a realization of a random 
process ^ such that F\I> = and 

* e w = : F f (RitUt), V (t)) < 1} (4) 

J to 

and 

/(■)€G :={/(•): f T (Q(t)f(t),f(t))<l}, (5) 

where Q(t) G l mxm , Q = Q' > 0, R(t) e M pxp , 
R' = R > and Q(t),R(t), R' 1 ^), Q~ l {t) are continuous 
functions of £ on [£o,F]. 

Suppose y(£) is observed in (3) for some x(-), / € G and 
^ G VF. The purpose of this section is to construct an 
algorithm with the following property: given a realization 
y(t), t € [to, T) of the random process Y = Hx(-) + ^, the 
algorithm produces an estimation of a linear function 

x(-)^(£,Fx(T)) 

having minimum mean-squared worst-case estimation er- 
ror. In what follows we will refer this algorithm as an a 
priori minimax mean-squared estimation in the direction I 
(^-estimation). Taking into account linearity of (3) we will 
look for ^-estimation among linear functions of y(-). Let us 
summarize the above discussion by rigorous mathematical 
definitions. 

Definition 2. Given u(-) e L 2 (t , T, W) and /el™ define 
a mean-squared worst-case estimation error 4 

a(T, u, £) := sup {E[(£, Fx(T)) - u{y)f : 

x{-),f {■)*{■) (6) 
(Fx) t - Gx + /, Fx(t ) = 0, /(•) g G, 77(0 G ^} 

3 8(G, /) = if / e G and +00 otherwise. 

4 Here the sup^.j means that we take the upper bound over all 

/■T 
{R(t)n(t),ri(t)) < 1 for all 

realizations rj(-) of 



2.1 Generalized Kalman Duality Principle 

The definition of the ^-estimation and error generalizes 
the notion of the linear minimax a priori mean-squared 
estimation, introduced in Nakonechny (1978). In order 
to find the ^-estimation we will follow a common way of 
deriving the minimax estimation Nakonechny (1978): first 
step is to obtain the expression for the worst-case error by 
means of the suitable duality concept, that is to formulate 
a dual control problem; next step is to solve it and to derive 
the minimax estimation. 

Next theorem generalizes the celebrated Kalman duality 
principle Brammcr and Siffling (1989) to non-causal DAEs. 

Theorem 3. (Generalized Kalman duality). The terror is 
finite iff 

(F'z) t (t) - -G'(t)z(t) + H'(t)u(t), F'z(T) = F'£ (7) 

(7) has a solution z(-). In this case the problem a(u) — > 
inf„ is equal to the following optimal control problem 

I(u) =mm{[ (Q^ 1 (z - v),z - v)dt\ 
v It 

t ( 8 ) 
+ / (R~ 1 u, u)dt — > min, 
Jto 

provided z(-) obeys (7) and v(-) obeys homogeneous (7). 

An obvious corollary of the Theorem 3 is an expression for 
the minimax observable subspace 

£(T) = {£ G K" : 3u(-),z(-)} (F' z) t +C' z-H'u = 0, F'z(T) 

Proof. Take £ G R" and u(-) G h 2 (t ,T,W) and suppose 
terror is finite. There exists some w(-) G L 2 (to,F, M m ) so 
that F'w(-) G H^to.^R") and F'w(T) = F'l A trivial 
example is w(t) = I. It was proved in Zhuk (2007) that 

(F'w{T),F+Fw{T)) - (Fx{t ),F+Fw(t )) = 

f T (9) 
/ ((Fx) u w) + ((F'w) u x)dt 

Jto 

if Fx(-) G H 1 (t ,T,R ,n ) and F'w(-) G H^to.T.R"). 
Noting that F = FF+F and using (9) and (3) one derives 

(£,Fx(T)) = {F'£,F+Fx(T)) = (F'w(T),F+Fx(T)) 
= [ ((Fx) t ,w) + {(F'w) u x)dt 

= f (f,w) + ((F'w) t + C'w,x)dt, 

Jto 

(10) 

Combining (10) with Erj = 0we have 



E[(l,Fx{T)) - u(y)f 



= [{£, Fx(T)) - f (H'u, x)dt} 2 + E[f (u{t),r]{t))dt] 
= [ f (/, w) + ((F'w) t + C'w - H'u, x)dt] 2 

J t a 

+ E[f T (u(t), V (t))dt} 2 

Jtn 



(11) 



Combining (4) with Cauchy inequality we obtain 



supE( (u,-q)dt) 2 = / (R- 1 u,u)dt (12) 
v J t J t 

(12) and a(u) < +oo imply the third line of (11) is 
bounded. Noting that J^(f,w)dt is bounded indepen- 
dently of x(-) one derives 

i-T 



sup{ f ((F'w) t +C'w-H'u, x)dt : (Fx) t = Cx+f, /(•) G G} 

x(-) Jt 

(13) 

It was proved in Zhuk (2009b) that 

sup {(C,x),LxeG}= inf {c(G, b), L'b = £} (14) 

x£3>(L) b£®(L>) 

provided 2>(L) := {x(-) G L 2 (i ,T,R n ) : Fx(-) G 
H^ic^R"),^^) =0} and 

(Lx){t) = (Fx) t (t) - C(t)x(t),x(-) € 9(h) (15) 

It was proved in Zhuk (2007) that 3>{L') := {b G 
h 2 (t ,T,R m ) : F'b(-) G H 1 ^, T, R m ), F'b(T) = 0} and 

L'b(t) = -(F'b) t - C'(t)b(t),b(-) € 9{U) (16) 

provided L is defined by (15). Setting C := (F'w)t + C'w — 
H'u we see from (13) that the right-hand part of (14) is 
finite. Using (16) one derives 

mi{c(G,b),-(F'b) t -C'(t)b(t) = (F'w) t +C'w-H'u} < +oo 

(17) 

Thus (17) implies 

(F'z) t + C'z = H'u, F'z = F'l 

with z := (w + b), &(•) G Sd(L'). This proves (7) has a 
solution z(-). Using integration-by-parts formula (9) and 
E'fy = and (12) one derives easily 

a(u) = sup / (w, f)dt+ / {R- 1 u,u)dt (18) 

/GGi J t J to 

with Gi denoting all /(•) e G such that (3) has a solution 
*(■)■ 

On the contrary, if z(-) is some solution of (7) then one 
derives (18) as it has been already done above. Therefore, 
there are only two cases: f-error is infinite or (18) holds. 

Note that 

sup(/ (f,z)dtf=sup{(f,z),feGnR(L)} 2 (19) 
fed Jto 

where R(L) is the range of the linear mapping L defined 
above by the rule (15). It was proved in Zhuk (2009b) that 

sup{(/,z),/ e G n R{L)} = M{c(G,z- v),v e N(L')} 

(20) 

provided intG n R{L) ^ 0. It is easy to see that the 
latter inclusion holds for L and G defined by (15) and (5) 
respectively. Recalling the definition of L' (formula (16)) 



and noting c 2 (G, z — v) = J ta (Q 1 (z — v), z — v)dt we derive 
: from (18)-(20) 

cr (u) = min{ f (Q^ 1 {z - v), z - v)dt} + j {R^ 1 u 1 u)dt 

v Jto Jto 

where L'v = 0. This concludes the proof. 



2.2 Optimality conditions and estimation algorithms 

Theorem 3 states that minimax estimation problem is 
equal to some optimal control problem for £ € C(T), which 
is called dual control problem. In the next proposition 
we introduce an approximate solution to the dual control 
problem without restricting the matrices F and C(t). 

Proposition 4- [Tikhonov rcgularization] Let I € £(T). 
For any e > the Euler-Lagrange system 

(F'z) t (t) = -C'(t)z(t) + H'{t)u + p, 
< 4&p)t(t) = C(t)p(t) + eQ-\t)z{t), (21) 

eu = Rp, Fp(t ) = 0, F'z(T) + F+Fp(T) = F'l 
has a unique solution p(e),z(e). and 

u(e) := -RHp(e) -> u in h 2 (t ,T, 
z{e) -> z in L 2 (t ,T,M m ), 



r(£) = lun -((F'£-F+Fp(T),Fp(T))- f \\p(e)\\ 2 dt), 

e Jto 

(22) 



where u, z denotes the solution of 

■ T 



min{min{ / (Q 1 (z — v),z — v)dt} + / (R 1 u,u)dt}, 

V Jto Jto 



(F'z) t (t) = -C'(t)z(t) + H'(t)u(t),F'z(T) = F'£ 



(23) 



Proof. Let terror be finite. Then (7) has a solution due 
to Theorem 3. Define (Hu) = (H'u,0), 1= (0,£) and set 
(Dz) = (-(F'z) t - C'z, F'z{T)) for z(-) e &{D) = {«(■) : 
F'z(-) e B 1 (t , T,R n )}. It is not difficult to see that the 
solution to (23) coincides with the solution (u, z) of the 
optimization problem 5 

||u|| 2 + ||z|| 2 ^min, D'z + Hu = l (*) 

This observation allows to apply the Tikhonov regular- 
ization Tikhonov and Arsenin (1977) method in order 
to derive (22). For simplicity assume that Q and R are 
equal to the identity mapping. Let us introduce Tikhonov 
function 

T £ {u,z) := \\F'z(T)-F'£\\ 2 + C \\(F'z) t + C'z — H'u\\ 2 dt 

Jto 

+ e f \\u\\ 2 + ||z|| 2 ^ = \\Dz - Hu l\\ 2 + e(\\u\\ 2 + \\z\\ 2 ) 

Jto 

(24) 

It is strictly convex and coercive. Thus it's minimum 
is attained at the unique point (u(s),z(s)). Moreover, 
(u(e),z(e)) goes to (u,z) in h 2 (to,T) as it follows from 
properties of the Tikhonov function Zhuk (2007). To 
conclude the proof it is sufficient to show that (u(e), z(e)) 



The norm is defined by (u, z) — > f T {Q 1 z,z) + (R 1 u,u)dt 



verifies (21). Using the argument of Zhuk (2007) we derive 
the Euler-Lagrange equation for (u(e), z(e)): 



Dz + Hu + p = I, 
D'p = ez, H'p = eu 



(25) 



where p = (p, q), H'p = [H'p, 0), D' is defined by the rule 
D'p = (Fp) t - Cp with p e 3>{D') = {p = (p, q) : Fp G 
H 1 (t , T, R m ),Fp(t ) = 0,q = F+ Fp{T) + d, Fd = 0}. For 
the detailed derivation of D' we refer the reader to Zhuk 
(2007). Now, introducing the definitions of D, D', H and 
H' into (25) we obtain (21). This proves the existence and 
uniqueness. Wc conclude with proving the last line in (22), 
which follows from (21) and the formula 

\\u(e)\\ 2 + \\z(eW^\\u\\ 2 + \\z\\ 2 =a(£) 

Remark 5. In fact, the above Proposition claims that £- 
estimation u is approximated by u{e) for any direction 
I € JC(T), provided u{e) is a linear transformation of a 
solution of Euler-Lagrange system (21) for the Tikhonov 
function (24) and (24) approximates the minimal worst- 
case error a(T, £). 

Now we will derive the suboptimal worst-case recursive 
estimator, acting 6 on a minimax observable subspacc. To 
do so we will introduce a splitting of (21) into differential 

and algebraic parts. 

Let D = diag(Ai . . . A r ) where Aj, i — 1, r := rangF are 

positive eigen values of FF' and set A = ( D 2 o rxn -r )_ 

Then Albert (1972) there exist S L € R mxm , S R € R nxn 
such that 

F = S L AS R , S L S' L = I, S R S' R = I, (26) 

Transforming (3) according to (26) and changing the 
variables one can reduce the general case to the case F = 
(I g). We split C(t), Q(t) and H'(t)R(t)H(t) according 
to the structure of F as follows: C{t) = (g* §), Q = 

(qI £)'™=(§ Define 
A(t) = (C' 3 Q^C 4 + S 2 ), B(t) = (C 2 - C' 4 Q^Q' 2 ), 
C(t, e) = -C[ + C' 3 Q^Q' 2 + A{t)M(t, e)B(t), 
W(t,s) = (eI + S± + C' i Q4 1 C 4 ),M(t,e) = W- 1 (t,e) 

Q(t, e) = ~A(t)M(t, e)A'(t) + 1 + + C' 3 Q?C 3 ], 
S(t, e) = eQi - eQ 2 Qi 1 Q' 2 + eB'(t)M(t, e)B{t), 
Corollary 6. [suboptimal estimation on a subspace] Let 

x t = -C'{e, t) - K{e, t)Q{e, t))x + -^H'Ry, 

z 1 =C(s,t)z 1 +Q(e,t)Kz 1 ,z 1 (T) = (I + K)~ 1 li (27) 
K t = -KC{e, t) - C'(e, t)K - KQ(e, t)K 
+ S{e,t),K(0) = 0,x(0) = 

with 7 I = (£ u l 2 ) and $(t, e ) = ( M{t , e)[£ K B % S l A ' { t )K] )- 
Then 



6 giving the estimation of the projection of the state vector x(t) onto 
a minimax observable subspace C(t) for all t £ [to,T]. 

7 £ is splitted into parts according to the splitting of R n induced by 
the block structure of F. 



sup E[(£, Fx{T)) - ((I + K(T, e))- 1 ^, x(T, e))] 2 
inf sup E[£(x) - u(y)} 2 , 

u xoJ,V 

&(£) = inf - [((/ + K(T, e))- 1 ^ K(T, e)(I + K(T, e))- 1 ^) 
e>0 e 



f \m,s) Zl \\ 2 dt\ 

Jtn 



(30) 



Proof. The idea 8 is to split the Euler-Lagrange sys- 
tem (21) into differential (pi,2i) and algebraic (p2, z 2) 
parts using the splittings of F, C, Q and H'RH, intro- 
duced above. We have 

Pi = C-tpi + C2P2 + e(Qizi + Q2Z2),Pi(to) = 0, 

z\ = -C[zi - C' 3 z 2 +pi + - {Stpi + S2P2), 

e 

® = C 3Pl + C A p 2 +e{Q' 2 z 1 + Q A z 2 ), (28) 

= -C 2 zi - C' 4 z 2 + - £ S' 2Pl + ^S 4 )P2 
z 1 (T)+p 1 (T)=£ 1 
Solving the algebraic equations for (p2,z 2 ) 

Z2 = Q-Ai-Q'2 - C 4 M)B Zl + ~(C 4 MA' ~ C 3 )pi], (0Q , 

p 2 = eMBz\ - MA' px, 

and substituting the resulting expressions to into differen- 
tial equations for (pi,z\) one obtains 

z 1 =C(e,t)z 1 + Q{e,t) Pl ,z 1 {T)+p 1 (T)=£ 1 , 

pi = -C'(e,t) Pl +S(e,t)z 1 , Pl (t a ) = 

Applying simple matrix manipulations one can prove that 
Q{e, t) > and S(e,t) > for e > 0, implying (30) is a 
non-negative Hamilton system for any e > 0. Therefore 
it is always solvable and the Riccati equation (27) has a 
unique symmetric non-negative solution. Note, that the 
unique solvability of (30) is also implied by Tikhonov 
method: (30) is equivalent to the Euler-Lagrange sys- 
tem (21), which is uniquely solvable. Now, by direct cal- 
culation we derive from (29)- (30) 

P= (Pi,P2) T = $(t,e)z 1 
Recalling that (Proposition 4) 

[ (y,-RHp(e))dt^ [ (y,u) in h 2 {t ,T,W),e -> 0, 

Jto e Jto 

and p = ®(t,e)z u h = C(e,t)z x + Q(e,t)Kz u z^T) = 
(I + K)~ 1 £i, we derive, integrating by parts, that 

[ T (y, -RHp(e))dt = f \ 1 <S>> H> 'Ry, Zl )dt = 

Jto £ Jto e 

((/ + i.)- 1 ^ 1 ,.i(r, £ )) 

where Xt is defined in (27). In the same manner we derive 
the expression for minimax error, recalling (22). This 
concludes the proof. 

Now we consider one special case when the DAE is regular 
and there is a possibility to derive the optimal state 
estimation algorithm. Let P 2 = P (V 2 = V) and R{P) = 
R(F) (R(V) = R(F% 



The same idea was used in Eremenko (1980). 



Proposition 7. [^-estimation and error] Let R((I—V)C'P) C 
R((I - V)C'(I - P)). Then for any !el" 

(Fp) t (f) - <7(t)p(t) + Q- 1 (t)z(t), Fp(t ) = 0, 
(F'z) t {t) = -C'(t)z(t)+H'(t)R(t)H(t)p(t),F'z(T) = F'l 

(31) 

has a solution. If p(-) and z(-) are some solution of (31) 
then, the ^-estimation u is given by u — RHp and the 
terror is represented by a{T,£) = (F'+£,Fp(T)). 

Proof. As above (26) we split C(t) and H' (t)R(t)H(t) 
and Q according to the structure of F. In this case (31) 
reads as 

pi = C\P\ + Q\Z\ + C 2 P2 + Q2Z2,Pl(to) = 0, 

z\ = -C[zi + Sipi - C' 3 z 2 + S2P2, z\{T) = ix, 
= C 3 pi + C iP2 + Q' 2 zi + Q±z 2 , 
= -C' 2Zl - C' 4 z 2 + S' 2 pi + S iP2 

Since Q" 1 > it follows that Q4 > implying z 2 = 
-Ql 1 (C 3 Pi + C 4 p 2 + Q' 2 zi) so that 

W(t,0)p 2 = B(t)z 1 -A'(t)p 1 (33) 

where A, B, W were defined above. It is easy to see that 
for our choice of F the proposition's assumption implies 
R(Cj) C R(C 2 ). Therefore (33) is always solvable (in the 
algebraic sense) and one solution has the form 

p 2 = W+(t,0)(B{t)z 1 -A'(t)p 1 ) 
Now we have to assume that p 2 € L 2 (f ,T). Substituting 
the representation for p 2 into (32) and noting that C±(I — 
W+(t, 0)W(t, 0)) = wc obtain 

pi = C+ (t)pi + S+ (t) Zl , Zl (T) = £ x 

zi = -C' + (t) Zl + Q+(t)pi,pi(t ) = (34) 

Pi(t o ) = 0,z 1 {T)=£ 1 

where C+{t) := d - Q 2 Q^C 3 - B'W+(t,0)A', S+(t) := 
Qi~Q 2 Q^Q 3 + B'W+(t,0)B, Q+{t) := Sx + C'gQ^Cz- 
AW + {t, Q)A' . Applying simple matrix manipulations one 
can prove that S + > 0, Q+ > so that (34) is a non- 
negative Hamilton system. Therefore it is always solvable. 

With help of (31) one easily shows I(u) — > and 
I(u) = a=(F'+£,Fp(T)). ' 

Remark 8. It is interesting to note that -S(t,e) — > <S+(f), 
sQ(t,e) -> Q+(t) and —C'(t,e) -> C + (t), provided e \ 
and the assumptions of the Proposition 7 hold. 
Corollary 9. [minimax estimation on a subspace] Let 
K = C+(t)K + KC+(t)' - KQ+{t)K + S+(t),K(t ) = 
where C + , Q + and S+ are defined above (proof of the 
Proposition 7). Then 

T {u lV ) = (£ u x{T)),&{T,£) = (K(T)£i,£i) 



J 

J to 



where x(to) = and 

x t = (C+(t)-KQ+(t))x + (K (B'-KA)w+(t,o))HRy(t) 

Proof. By direct calculation one finds that Kzi , Z\ ver- 
ify (34) so that pi = Kz\. Combining this and (33) with 
u = RHp and a(T,£) = (F'+£, Fp(T)) (Proposition 7) one 
obtains the statement of the corollary. 

2.3 Numerical example: non-causal non- stationary DAE 



Let 



F=(l ° ),C(t) = ( - { \ ) l),H(t) = (o 1) 

Then dct(F - XC(t)) = if c 3 (t) = 0. The corresponding 
DAE reads 

x\ = -xi + x 2 + fi(t), 
= c 3 (f)x 1 (f) + / 2 (f),z 1 (0) = 

Set /1 = for simplicity. We have x\(t) = J * cxp(s — 
t)x 2 (s)ds and 

Cs(t) / exp(s - t)x 2 (s)ds = -f 2 (t) (36) 



(35) 



Set cj(t) = if c 3 (f) = and 
formally 



otherwise. Then, 



ar 2 (t) = exp(-f)-(-exp(f)c+(f)/ 2 (f)) + <;(f), (37) 

with c 3 (i) j* * exp(f - s)w(s)d.s = for f e [0,T]. We see 

that f 2 must be able to suppress the growth of c 3 near 
points where 03(f) = 0, in order to (37) belong to L 2 (0, T). 
Taking f 2 = exp(-c 3 f )b, &(■) G H^fo,! 1 ) we obtain 
/2(f) € K{c 3 {t)) and c+/ 2 € H^O.T). Therefore = 
J exp(s— t)x 2 (s)ds andx 2 (f) = — c^f) exp(— C3 (f))(6(f)+ 

^ (*)) - (4W ex P( -c 3~(*)))« & (*) + From this formula 
we see that x 2 is driven by v only if 03(f) = 0. If c 3 (f) 7^ 
then x 2 is driven by / 2 , its derivative and w. From the 
analytical point of view this implies the corresponding 
DAE is ill-posed: x\ is non-unique and is not continuous 
with respect to input data. Namely, as the differential 
operator is unbounded in L 2 (fo,T), we see that x 2 is not 
continuous with respect to / 2 , implying ill-posedness. As 
x 2 depends on an arbitrary function v from some linear 
subspace we have non- uniqueness. 

The aim is to estimate x\(t), provided y(t) = x 2 (t) + 77(f) 
is measured and 

E - V 2 (t)dt<l, Er?(-)=0 (38) 

and (x\,x 2 ) obeys 

xi = -x 1 +x 2 + fx, a;i(0)=0 
= c 3 (t) Xl (t) + f 2 (t) 

with f 2 (f) = exp(-c+(f))6(f), &(•) e H^O.T) and 
Jo 



< > eX ^f 2 dt < 1 



(40) 



Figure 1 describes the observations y of x 2 , perturbed 
by the non-Gaussian noise r/(-), arbitrary function v(-) 
and uncertain / 2 (-), provided rj(-) verifies (38), / 2 (-) ver- 
ifies (40) and c 3 (f) J * exp(f - s)v(s)ds = for f e [0,T]. 
The sub-optimal estimation x(f,e) and error are given by 

x t (f,e) =[ - 1 - K(l + -M^ )]x(f,e) + ^L,x(e,0) = 0. 



eexp(v / f) ' 



6 + TV 



Jf t = -2K - (1 + ^r)i^ 2 + e(l + -g ), ^(0, e) = 0, 



eexp(v / f) 
A-(t,e) 



£ ) + 76 



+e) 2 



:= e[(X(f, £) + 

z t £ (f;r) = [l + if(l + - 

eexp(vf) 

Note that the pencil sF' - C" - iJ'i? is regular. As (31) 
is solvable in this case we apply Proposition 7 in order to 



}z^s;t) 2 ds], 



2C ' W )}z^t;r),z%T,T) = (K(T,e) + l)- 1 




y(t) 
* 2 (t) 



Fig. 1. Observation's noise 








— m 




— i(t,e) 


: 


-- o-(t,e) 




o-(t) 







Fig. 2. Optimal estimation and error, suboptimal estima- 
tion and error and x\ 

derive the minimax estimation x and minimax estimation 
error <r(i): 

i t =[-l-K(t)—f±-)]x + y(t), 5(0) = 0, 
cxp(Vt) 

K t = -2K - —fjfr)K 2 + (1 + -), K(0, s) = 0, 
exp(Vt) 6 

a 2 (Mi) :=A-(t)^ 
In Figure 2 the comparison of the optimal estimator x and 
error a(t) with sub-optimal estimator x(e) and error <r(i, s) 
are presented, provided £ = exp(— 100) and l\ = 0. 
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